Back

serac: a R package for ShortlivEd RAdionuclide Chronology of recent sediment cores

Rosalie Bruel
April 5th, 2019

Content


I had to produce some age-depth model plot for my PhD thesis. My colleagues working on dating would usually use a spreadsheet, and then a paying photo editor to edit their figures. I do not have the license for either of these, so of course I turned to R. My colleague Pierre Sabatier thought it was great and that we should turn the code in something even more automatic.

1. Get your data

We are sharing within the package data from Lake Allos (France), so users can follow along. The publication by Wilhelm et al (2012) is accessible here.

Lake Allos (source: Flickr)left: Lake Allos (Source: Flickr); right: Location, hydrological, geological and geomorphological characteristics of Lake Allos and its catchment (Fig. 1 in Wilhelm et al 2012)

Your dataset should follow a similar format than the one below.

2. Download serac

You can download serac from my GitHub repository using the packages devtools.

3. Test your hypotheses

At the minimum, your function will be:

serac(name="ALO09P12",coring_yr=2009)

 Sedimentation rate (CFCS model): V= 1.197mm/yr, R2= 0.9597
                          Error:     +/- 0.053mm/yr
 Inventory (Lead): 11707.91 Bq/m2 (range: 11176-12240 Bq/m2)
 Inventory (Cesium): 21549.328 Bq/m2 (range: 21405-21694 Bq/m2)



 ________________________

 The calculation took 0.856 seconds.

Several sedimentation hypotheses were tested, and this is what the author chose as the best model:

serac(name="ALO09P12",coring_yr=2009,model=c("CFCS"),plotphoto=FALSE,minphoto=c(0),maxphoto=c(210),plot_Pb=T,plot_Am=T,plot_Cs=T,Cher=c(30,40),Hemisphere=c("NH"),NWT=c(51,61),sedchange=c(75.5),plot_Pb_inst_deposit=T,inst_deposit=c(20,28,100,107,135,142,158,186),suppdescriptor=TRUE,descriptor_lab=c("Ca/Fe"),historic_d=c(20,28,100,107,135,142,158,186),historic_a=c(1994,1920,1886,1868),historic_n=c("sept1 994 flood","1920 flood","1886 flood","1868 flood ?"), min_yr=c(1750),dmax=c(180), plotpdf=TRUE,preview=F)

 Sedimentation rate (CFCS model) 0-75.5mm: V= 1.096mm/yr, R2= 0.9451
                          Error:     +/- 0.088mm/yr

 Sedimentation rate (CFCS model) 75.5mm-bottom: V= 0.796mm/yr, R2= 0.9226
                          Error:     +/- 0.082mm/yr
 Approximation of age at change(s) in sedimentation rate:
     Best Age: 1947 (incertitude: 1942-1952)


 Age approximation of instantaneous deposit(s) from CFCS model:
     The instantaneous deposit at 20-28 mm has an estimated range of: 1989-1992.

     The instantaneous deposit at 100-107 mm has an estimated range of: 1905-1928.

     The instantaneous deposit at 135-142 mm has an estimated range of: 1866-1897.

     The instantaneous deposit at 158-186 mm has an estimated range of: 1844-1879.

 Inventory (Lead): 10287.66 Bq/m2 (range: 9828-10748 Bq/m2)
 Inventory (Cesium): 19023.688 Bq/m2 (range: 18895-19153 Bq/m2)


 ________________________

 The calculation took 0.443 seconds.

Source

Table of content created using this method.

is loading comments…
LS0tCm91dHB1dDoKICBodG1sX2RvY3VtZW50OgogICAgZGZfcHJpbnQ6IHBhZ2VkCi0tLQoKPCEtLSBHbG9iYWwgc2l0ZSB0YWcgKGd0YWcuanMpIC0gR29vZ2xlIEFuYWx5dGljcyAtLT4KPHNjcmlwdCBhc3luYyBzcmM9Imh0dHBzOi8vd3d3Lmdvb2dsZXRhZ21hbmFnZXIuY29tL2d0YWcvanM/aWQ9VUEtMTM0ODcwODA4LTEiPjwvc2NyaXB0Pgo8c2NyaXB0PgogIHdpbmRvdy5kYXRhTGF5ZXIgPSB3aW5kb3cuZGF0YUxheWVyIHx8IFtdOwogIGZ1bmN0aW9uIGd0YWcoKXtkYXRhTGF5ZXIucHVzaChhcmd1bWVudHMpO30KICBndGFnKCdqcycsIG5ldyBEYXRlKCkpOwoKICBndGFnKCdjb25maWcnLCAnVUEtMTM0ODcwODA4LTEnKTsKPC9zY3JpcHQ+CgpgYGB7ciBpbmNsdWRlPUZBTFNFfQpzb3VyY2UoInJlbmRlcl90b2MuUiIpCmBgYAoKCjxzdHlsZT4KLmZ1bGwtcGFnZXsKICBmbG9hdDogbGVmdDsKICB3aWR0aDogNzAlOwogIHZlcnRpY2FsLWFsaWduOiB0b3A7CiAgdGV4dC1hbGlnbjogbGVmdDsKICBtYXJnaW4tbGVmdDogMTUlOwogIG1hcmdpbi1yaWdodDogMTUlOwp9Ci5ib2R5IHsKICBiYWNrZ3JvdW5kLWltYWdlOiB1cmwoIkNvcmVzL0FMTzA5UDEyL0FMTzA5UDEyLmpwZyIpOwogIGJhY2tncm91bmQtcmVwZWF0OiBuby1yZXBlYXQ7CiAgYmFja2dyb3VuZC1hdHRhY2htZW50OiBmaXhlZDsKICBiYWNrZ3JvdW5kLXNpemU6IGNvbnRhaW47CiAgYmFja2dyb3VuZC1wb3NpdGlvbjogMTUlOwp9Ci5jb250ZW50IHsKICAgIHBhZGRpbmc6MjBweDsKICAgIHdpZHRoOjcxMHB4OwogICAgcG9zaXRpb246cmVsYXRpdmU7CiAgICBiYWNrZ3JvdW5kOiByZ2IoMjA0LCAyMDQsIDIwNCk7IC8qIEZhbGxiYWNrIGZvciBvbGRlciBicm93c2VycyB3aXRob3V0IFJHQkEtc3VwcG9ydCAqLwogICAgYmFja2dyb3VuZDogcmdiYSgyNTUsMjU1LDI1NSwgMC44NSk7Cn0KPC9zdHlsZT4KCjxhIGhyZWY9Ii4uL2Jsb2cuaHRtbCI+IEJhY2sgPC9hPgoKPGRpdiBjbGFzcz0iYm9keSI+Cgo8ZGl2IGNsYXNzPSJmdWxsLXBhZ2UgY29udGVudCI+CgojIyBzZXJhYzogYSBSIHBhY2thZ2UgZm9yIFNob3J0bGl2RWQgUkFkaW9udWNsaWRlIENocm9ub2xvZ3kgb2YgcmVjZW50IHNlZGltZW50IGNvcmVzCjxkaXYgc3R5bGU9InRleHQtYWxpZ246cmlnaHQiPgpSb3NhbGllIEJydWVsICAKKkFwcmlsIDVedGheLCAyMDE5Kgo8L2Rpdj4KCl9fQ29udGVudF9fCmBgYHtyIHRvYywgZWNobz1GQUxTRX0gCnJlbmRlcl90b2MoIlR1dG8yX3NlcmFjLlJtZCIpCmBgYApgYGB7ciBnZXQgZGF0YSBjYW1wLCBtZXNzYWdlPUZBTFNFLCB3YXJuaW5nPUZBTFNFLCBpbmNsdWRlPUZBTFNFLCBwYWdlZC5wcmludD1GQUxTRX0KaWYoIXJlcXVpcmUoZGV2dG9vbHMpKSB7CiAgaW5zdGFsbC5wYWNrYWdlcygiZGV2dG9vbHMiKQp9CmRldnRvb2xzOjppbnN0YWxsX2dpdGh1YigiZGF0YWNhbXAvdHV0b3JpYWwiKQpgYGAKYGBge3IsIGluY2x1ZGU9RkFMU0V9CnR1dG9yaWFsOjpnb19pbnRlcmFjdGl2ZSgpCmBgYAogIAoqKiogIAoKCkkgaGFkIHRvIHByb2R1Y2Ugc29tZSBhZ2UtZGVwdGggbW9kZWwgcGxvdCBmb3IgbXkgUGhEIHRoZXNpcy4gTXkgY29sbGVhZ3VlcyB3b3JraW5nIG9uIGRhdGluZyB3b3VsZCB1c3VhbGx5IHVzZSBhIHNwcmVhZHNoZWV0LCBhbmQgdGhlbiBhIHBheWluZyBwaG90byBlZGl0b3IgdG8gZWRpdCB0aGVpciBmaWd1cmVzLiBJIGRvIG5vdCBoYXZlIHRoZSBsaWNlbnNlIGZvciBlaXRoZXIgb2YgdGhlc2UsIHNvIG9mIGNvdXJzZSBJIHR1cm5lZCB0byBSLgpNeSBjb2xsZWFndWUgUGllcnJlIFNhYmF0aWVyIHRob3VnaHQgaXQgd2FzIGdyZWF0IGFuZCB0aGF0IHdlIHNob3VsZCB0dXJuIHRoZSBjb2RlIGluIHNvbWV0aGluZyBldmVuIG1vcmUgYXV0b21hdGljLgoKIyMjIyAxLiBHZXQgeW91ciBkYXRhCldlIGFyZSBzaGFyaW5nIHdpdGhpbiB0aGUgcGFja2FnZSBkYXRhIGZyb20gTGFrZSBBbGxvcyAoRnJhbmNlKSwgc28gdXNlcnMgY2FuIGZvbGxvdyBhbG9uZy4gVGhlIHB1YmxpY2F0aW9uIGJ5IFdpbGhlbG0gKmV0IGFsKiAoMjAxMikgaXMgYWNjZXNzaWJsZSA8YSBocmVmPWh0dHBzOi8vd3d3LnNjaWVuY2VkaXJlY3QuY29tL3NjaWVuY2UvYXJ0aWNsZS9waWkvUzAwMzM1ODk0MTIwMDAyOTQgdGFyZ2V0PSJfYmxhbmsiPmhlcmU8L2E+LiAKCjxkaXYgc3R5bGU9IndpZHRoOjEwMCU7IG1hcmdpbi1yaWdodDogNXB4OyBtYXJnaW4tYm90dG9tOiA1cHg7IG1hcmdpbi10b3A6IDVweDsgZm9udC1zaXplOjgwJTsgdGV4dC1hbGlnbjpjZW50ZXI7Ij48aW1nIHNyYz0iaW1hZ2VzL1R1dG8yLzFfQWxsb3NfRmxpY2tyX1dpbGhlbG0yMDEyLnBuZyIgYWx0PSJMYWtlIEFsbG9zIChzb3VyY2U6IEZsaWNrcikiICBzdHlsZT0icGFkZGluZy1ib3R0b206MC41ZW07IiBhbGlnbjoicmlnaHQiICAvPl9fbGVmdF9fOiBMYWtlIEFsbG9zIChTb3VyY2U6IEZsaWNrcik7IF9fcmlnaHRfXzogTG9jYXRpb24sIGh5ZHJvbG9naWNhbCwgZ2VvbG9naWNhbCBhbmQgZ2VvbW9ycGhvbG9naWNhbCBjaGFyYWN0ZXJpc3RpY3Mgb2YgTGFrZSBBbGxvcyBhbmQgaXRzIGNhdGNobWVudCAoRmlnLiAxIGluIFdpbGhlbG0gZXQgYWwgMjAxMik8L2Rpdj4KCgpZb3VyIGRhdGFzZXQgc2hvdWxkIGZvbGxvdyBhIHNpbWlsYXIgZm9ybWF0IHRoYW4gdGhlIG9uZSBiZWxvdy4gCgoKCgpgYGB7ciByZWFkIGFuZCB2aXN1YWxpemUgZGF0YSwgZWNobz1GQUxTRX0KYWxvIDwtIHJlYWQudGFibGUocGFzdGUwKGdldHdkKCksIi9Db3Jlcy9BTE8wOVAxMi9BTE8wOVAxMi50eHQiKSwgaGVhZGVyPVQpCmFsbwpgYGAKCgoKIyMjIyAyLiBEb3dubG9hZCBzZXJhYwpZb3UgY2FuIGRvd25sb2FkICpzZXJhYyogZnJvbSBteSBHaXRIdWIgcmVwb3NpdG9yeSB1c2luZyB0aGUgcGFja2FnZXMgZGV2dG9vbHMuCgpgYGB7ciBsb2FkIHBhY2thZ2VzLCBldmFsPUZBTFNFLCBpbmNsdWRlPUZBTFNFfQojIGluc3RhbGwgc2VyYWMgd2l0aCB0aGUgcGFja2FnZSBkZXZ0b29scyAoZG93bmxvYWQgaXQgaWYgbm90IGFscmVhZHkgaW4gbGlicmFyeSkKIyByZXF1aXJlKGRldnRvb2xzKSAjIElmIGRldnRvb2xzIG5lZWQgdG8gYmUgZG93bG9hZGVkCiNsaWJyYXJ5KGRldnRvb2xzKQojZGV2dG9vbHM6Omluc3RhbGxfZ2l0aHViKCJyb3NhbGllYi9zZXJhYyIpCiMgbG9hZCB0aGUgc2VyYWMgbGlicmFyeQojbGlicmFyeShzZXJhYykKCiMgSGVscCBmdW5jdGlvbiBzdW1tYXJpemluZyB0aGUgZmlyc3Qgc3RlcHMKI2hlbHBfc2VyYWMoKQoKIyBVc2VyIGluZm9ybWF0aW9ucyAjIyMjCiMgRW50ZXIgc29tZSBpbmZvcm1hdGlvbiBvbiB0aGUgdXNlciAobmFtZSwgT1JDSUQsIGFmZmlsaWF0aW9uLCBlbWFpbCkuCiMgRG8gaXQgb25jZSwgYW5kIHRoZXNlIGluZm9ybWF0aW9uIHdpbGwgYXV0b21hdGljYWxseSBiZSBpbmNsdWRlZCB0byBhbnkgbmV3IHByb2plY3QuCiMgUnVuIHVzZXJfaW5mb3MoKSBpbiB0aGUgY29uc29sZSBvciBzaXAgdGhpcyBzdGVwCiMgdXNlcl9pbmZvcygpCiNnZXR3ZCgpCmBgYAoKCgoKIyMjIyAzLiBUZXN0IHlvdXIgaHlwb3RoZXNlcwpBdCB0aGUgbWluaW11bSwgeW91ciBmdW5jdGlvbiB3aWxsIGJlOgpgYGB7ciBydW4gYmFzaWMgbW9kZWwsIG1lc3NhZ2U9RkFMU0UsIHdhcm5pbmc9RkFMU0UsIHBhZ2VkLnByaW50PUZBTFNFfQojc2VyYWMobmFtZT0iQUxPMDlQMTIiLGNvcmluZ195cj0yMDA5KQpgYGAKClNldmVyYWwgc2VkaW1lbnRhdGlvbiBoeXBvdGhlc2VzIHdlcmUgdGVzdGVkLCBhbmQgdGhpcyBpcyB3aGF0IHRoZSBhdXRob3IgY2hvc2UgYXMgdGhlIGJlc3QgbW9kZWw6CmBgYHtyfQojc2VyYWMobmFtZT0iQUxPMDlQMTIiLGNvcmluZ195cj0yMDA5LG1vZGVsPWMoIkNGQ1MiKSxwbG90cGhvdG89RkFMU0UsbWlucGhvdG89YygwKSxtYXhwaG90bz1jKDIxMCkscGxvdF9QYj1ULHBsb3RfQW09VCxwbG90X0NzPVQsQ2hlcj1jKDMwLDQwKSxIZW1pc3BoZXJlPWMoIk5IIiksTldUPWMoNTEsNjEpLHNlZGNoYW5nZT1jKDc1LjUpLHBsb3RfUGJfaW5zdF9kZXBvc2l0PVQsaW5zdF9kZXBvc2l0PWMoMjAsMjgsMTAwLDEwNywxMzUsMTQyLDE1OCwxODYpLHN1cHBkZXNjcmlwdG9yPVRSVUUsZGVzY3JpcHRvcl9sYWI9YygiQ2EvRmUiKSxoaXN0b3JpY19kPWMoMjAsMjgsMTAwLDEwNywxMzUsMTQyLDE1OCwxODYpLGhpc3RvcmljX2E9YygxOTk0LDE5MjAsMTg4NiwxODY4KSxoaXN0b3JpY19uPWMoInNlcHQxIDk5NCBmbG9vZCIsIjE5MjAgZmxvb2QiLCIxODg2IGZsb29kIiwiMTg2OCBmbG9vZCA/IiksIG1pbl95cj1jKDE3NTApLGRtYXg9YygxODApLCBwbG90cGRmPVRSVUUscHJldmlldz1GKQpgYGAKCgo8aW1nIHNyYz0iQ29yZXMvQUxPMDlQMTIvQUxPMDlQMTIucGRmIiBzdHlsZT0id2lkdGg6MTAwJTsgbWFyZ2luLXJpZ2h0OiA1cHg7IG1hcmdpbi1ib3R0b206IDVweDsgbWFyZ2luLXRvcDogNXB4IiBhbGlnbj0iY2VudGVyIj4KCiMjIyMgU291cmNlICAKClRhYmxlIG9mIGNvbnRlbnQgY3JlYXRlZCB1c2luZyA8YSBocmVmPWh0dHBzOi8vd3d3LmdhcnJpY2thZGVuYnVpZS5jb20vYmxvZy9hZGQtYS1nZW5lcmF0ZWQtdGFibGUtb2YtY29udGVudHMtYW55d2hlcmUtaW4tcm1hcmtkb3duLyB0YXJnZXQ9Il9ibGFuayI+dGhpcyBtZXRob2Q8L2E+LgoKCjwvZGl2PgoKCgo8IS0tIGJlZ2luIHd3d3cuaHRtbGNvbW1lbnRib3guY29tIC0tPgogPGRpdiBpZD0iSENCX2NvbW1lbnRfYm94Ij48YSBocmVmPSJodHRwOi8vd3d3Lmh0bWxjb21tZW50Ym94LmNvbSI+PC9hPiBpcyBsb2FkaW5nIGNvbW1lbnRzLi4uPC9kaXY+CiA8bGluayByZWw9InN0eWxlc2hlZXQiIHR5cGU9InRleHQvY3NzIiBocmVmPSJodHRwczovL3d3dy5odG1sY29tbWVudGJveC5jb20vc3RhdGljL3NraW5zL2Jvb3RzdHJhcC90d2l0dGVyLWJvb3RzdHJhcC5jc3M/dj0wIiAvPgogPHNjcmlwdCB0eXBlPSJ0ZXh0L2phdmFzY3JpcHQiIGlkPSJoY2IiPiAvKjwhLS0qLyBpZighd2luZG93LmhjYl91c2VyKXtoY2JfdXNlcj17fTt9IChmdW5jdGlvbigpe3ZhciBzPWRvY3VtZW50LmNyZWF0ZUVsZW1lbnQoInNjcmlwdCIpLCBsPWhjYl91c2VyLlBBR0UgfHwgKCIiK3dpbmRvdy5sb2NhdGlvbikucmVwbGFjZSgvJy9nLCIlMjciKSwgaD0iaHR0cHM6Ly93d3cuaHRtbGNvbW1lbnRib3guY29tIjtzLnNldEF0dHJpYnV0ZSgidHlwZSIsInRleHQvamF2YXNjcmlwdCIpO3Muc2V0QXR0cmlidXRlKCJzcmMiLCBoKyIvanJlYWQ/cGFnZT0iK2VuY29kZVVSSUNvbXBvbmVudChsKS5yZXBsYWNlKCIrIiwiJTJCIikrIiZvcHRzPTE2ODYyJm51bT0xMCZ0cz0xNTUzMjczMjYyNDMyIik7aWYgKHR5cGVvZiBzIT0idW5kZWZpbmVkIikgZG9jdW1lbnQuZ2V0RWxlbWVudHNCeVRhZ05hbWUoImhlYWQiKVswXS5hcHBlbmRDaGlsZChzKTt9KSgpOyAvKi0tPiovIDwvc2NyaXB0Pgo8IS0tIGVuZCB3d3cuaHRtbGNvbW1lbnRib3guY29tIC0tPgoK